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We study two intertwined globally coupled networks of noisy Kuramoto phase oscillators that have 
the same natural frequency, but differ in their perception of the mean field and their contribution to 
it. Such a give-and-take mechanism is given by asymmetric in- and out-coupling strengths which can 
be both positive and negative. We uncover in this minimal network of networks intriguing patterns 
of discordance, where the ensemble splits into two clusters separated by a constant phase lag. If 
it differs from tt, then traveling wave solutions emerge. We observe a second route to traveling 
waves via traditional one-cluster states. Bistability is found between the various collective states. 
Analytical results and bifurcation diagrams are derived with a reduced system. 

PACS numbers: 05.40.-a, 05.45.Xt, 87.10.Ca 


I. INTRODUCTION II. MODEL 


From neurons in the brain to cells in the heart, from 
electrons in superconductors to planets in the universe, 
collective oscillations are ubiquitous [T]. Clearly, the un¬ 
derlying mutual synchronization critically relies on the 
presence of interactions. If individuals are allowed to 
interact non-uniformly, the collective behavior remains 
particularly elusive. Beyond incoherence and synchro¬ 
nization, what other collective states can emerge, and 
under which conditions? Motivated by Daido’s seminal 
work on “oscillator glasses” [5], Hong and Strogatz tack¬ 
led this question recently in a series of papers mm- They 
investigated two mutually globally coupled populations 
of Kuramoto phase oscillators that differed in their cou¬ 
pling strengths. Two scenarios of mixed attractive and 
repulsive interactions were distinguished. In the first case 
some oscillators’ phases repel the phases of all the oth¬ 
ers, while the remaining attract all the other phases [3]. 
This situation resembles neural networks with excitatory 
and inhibitory connections [S]. In the second case some 
oscillators tend to align with the mean field, while others 
oppose it favoring an antiphase alignment [1]. This ver¬ 
sion is analogous to sociophysical models of opinion for¬ 
mation [B]. Surprisingly, only the second scenario led to 
enriched dynamics beyond the traditional order-disorder 
transition. 

In this paper, we unify both coupling scenarios. We un¬ 
veil the existence of two routes to traveling waves, which 
in parameter space are surrounded either by states of 
diametral two-cluster synchronization or by one-cluster 
partially synchronous states. Differently than in Refs. 
mm we consider temporal fluctuations acting on the 
frequencies. Our main contribution is the derivation of 
bifurcation diagrams for all possible collective states. 


One of the most prominent models describing phe¬ 
nomena of mutual synchronization is due to Kuramoto 
[7]. It describes how the phases of coupled oscillators 
evolve in time, and is applicable to systems of nearly 
identical, weakly coupled limit-cycle oscillators. We con¬ 
sider a stochastic version with twofold disordered cou¬ 
pling strengths: 

K , 

^,{t) = Wo + ^ X ~ 

i=i 

where i = 1,..., A^. Using the notion of give-and-take 
as a metaphor, oscillator i contributes to the mean field 
with weight Gi and at the same time incorporates the 
mean activity with weight Ki into its own dynamics. Ac¬ 
cordingly, we call Ki in- and Gi out-coupling strength, 
respectively. Grouping together oscillators with the same 
coupling strengths, the number of different pairs {Ki, Gi) 
coincides with the number of subpopulations. 

We consider two equally-sized subpopulations denoted 
by “1” and “2.” Hence, the oscillators are distinguished 
by a pair of coupling strengths, (ATi, Gi) or {K 2 , G 2 ), and 
all of those can be positive or negative. We choose the 
parametrization 

Ari,2 = Ko±^—, Gi,2 = Go±^-. (2) 

Kq and Go are average in- and out-coupling strengths, 
while ART and AG give corresponding mismatches. If 
\AK\ /2 > iRTol or |AG| /2 > |Go|, then half of the cou¬ 
plings are positive (attractive) and half negative (repul¬ 
sive). In such cases we speak of mixed interactions. Note 
that leads to point symmetries, because changing 
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{Ko,Go) ^ {-Ko,-Go) or {AK,AG) ^ {-AK,-AG) 
yields the same situations. In the following, all oscillators 
have the same constant natural frequency Wq- There¬ 
fore, by virtue of the rotational symmetry, we can set 
Wo = 0 without loss of generality. Time-dependent dis¬ 
order ^i{t) is included as Gaussian white noise, = 

0, = 2DSijS{t — f). The angular brackets de¬ 

note averages over different realizations of the noise and 
the single non-negative parameter D denotes the noise 
intensity. The noise terms ^i{t) can be regarded as an 
aggregation of various stochastic processes [5]. In Refs. 
mm it is found for the deterministic case D — 0 that 
mixed out-couplings alone do not enable more than par¬ 
tial synchronization, whereas mixed in-couplings yield 
traveling waves reached through diametrically synchro¬ 
nized states. We intermingle both types of mixing and 
explore in particular whether traveling waves persist in 
the presence of noise D > 0. Without loss of generality 
all subsequent results are obtained with D = 0.5, but for 
illustration we keep D in the derivations. 

“Discordant synchronization” is used here as an um¬ 
brella term for situations where the ensemble splits into 
two partially synchronized clusters. This will include 
traveling waves and 7r-states, the latter being the extreme 
form of discordance with two oscillator populations anti¬ 
aligned to each other. 


III. THEORY 


We investigate the thermodynamic limit TV —> cx), 
where propagation of molecular chaos [Ij allows us to 
describe each population by a one-oscillator probabil¬ 
ity density ^' 1 ^ 2 )- Normalization 

requires = 1 V t. For given coupling 

strengths Ki ^2 and Gi, 2 , pip{4’,t)d4> denotes the fraction 
of oscillators with phase between </> and (j) -\- dcj) at time 
t. The densities are governed by the nonlinear Fokker- 
Planck equations m- 


dpi ,2 _ „ 5^1,2 
dt d(j)^ 


A 

d(j) 


[Ki,2Rs\n (0 - (j)) pi,2] 


( 3 ) 
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FIG. 1. (Color online) Snapshots of oscillators’ phases, (a) 
Stationary rr-state at Kq = 1, mean phase difference 5 = tt, 
no net flux; (b) Traveling wave at Kq = 3, flow downwards 
with stationary profile, indicated by black arrows; red dots 
depict phase values, blue lines of different lengths and direc¬ 
tions indicate instantaneous frequencies. Remaining parame¬ 
ters: AK = 3, Go = 2, AG = 10, D = 0.5. For visualization 
the two subpopulations are separated into two halves. 


r{t) = i \ri{t) _|_ r 2 {t) , which differs from 

Eq. Q by taking out Gi_ 2 - The order parameter values 
lie between zero (incoherence) and 1 (complete synchro¬ 
nization). The variables 0(t) and 0i,2(t) stand for the 
corresponding mean phases. Of special interest is the 
phase lag 5{t), i.e. the difference in the mean phases of 
the two populations, 5{t) = 0i(t) — 02(t)- 

In the sequel, variables without a dot or other indicated 
time dependence refer to the long-time limit. Let us out¬ 
line the four qualitatively different self-organized states 
observed here after some transient dynamics, (i) In the 
incoherent state the whole population of oscillators ro¬ 
tates asynchronously, ri ^2 = 0. (ii) The classical partially 
synchronized state has zero phase lag, ri _2 >0, (5 = 0. 
We use “zero-lag sync” as a shortcut to denote this state, 
(iii) The 7r-state describes a partially synchronized state, 
where the two subpopulations are anti-aligned to each 
other, ri _2 >0, 5 = tt. (iv) In the traveling wave state 
the whole population is also partially synchronized, but 
oscillates with a frequency different from the frequency 
of single oscillators. This spontaneous change in rhythm 
is induced by a phase lag that is neither zero nor tt, 
f’ 1,2 >0, 0 < (5 < TT. We calculate the wave speed as 


The global mean-field amplitude R{t) and phase 0(f) 
follow from a superposition: 


R(t)e 


,»©(*) — 


1 r 

2 


ri(f) Gi -br 2 (f) G 2 


( 4 ) 

Note that subpopulations of different sizes can be treated 
simply by rescaling Gi, 2 . The local mean-field variables 
obey 


p2'K 

2 ((/>', t). (5) 

^0 

The level of synchrony in the two subpopulations is mea¬ 
sured separately by ri, 2 {t), whereas for the global mea¬ 
sure we take the classical Kuramoto order parameter 


1 

^ ( 6 ) 

1=1 

where (...)* represents a long-time average [3]. 

In Fig. example snapshots from simulations of 
N = 10^ oscillators are shown {N = 10^ are equidis- 
tantly chosen for visualization). Figure a) displays a 
TT-state with ri = 0.98, r 2 = 0.89, r = 0.05, |( 5 | = tt and 
|n| = 0. Figure [^b) shows a traveling wave state with 
ri = 0.98, r 2 = 0.75, r = 0.80, \S\ = 1.42 and \n\ = 3.17. 
It is equally possible that the wave runs in the other di¬ 
rection, depending on initial phases and realization of the 
noise. Note that perfect synchrony, ri _2 = !> cannot be 
achieved with finite coupling strengths, if an infinitesimal 
amount of noise is present. 
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In order to analytically investigate the collective dy¬ 
namics that are governed by (§-(§, we approximate 
the phase distributions in the two populations by time- 
dependent Gaussians. This well-known method is mo¬ 
tivated by numerical observations HU HI]. It has been 
generalized recently to encompass heterogeneities in cou¬ 
plings unis]. Extending those derivations to the present 
case, we obtain the following three-dimensional system of 

ode’s: 

ri = - riD + [riGi -|- r 2 G 2 cos S] , 

r 2 = - r 2 D + [?’ 2 G 2 -I- riGi cos 5] , 

S = - -k rl) Kir2G2 + -k r^) K2riGi] . 

( 7 ) 

All the four aforementioned collective states are fixed 
points of Q with ri ^2 = <5 = 0. Two types of fixed 
point solutions have to be distinguished, because there 
are two possibilities that <5 = 0 holds: 

5 = mTT, TO G Z, (8) 

0 = (rf ^ -k rl) Kir 2 G 2 + -k r|) K 2 riGi. (9) 

Equation describes zero-lag and 7r-states, whereas 
Eq. § underlies traveling waves. Intermediate phase 
lags 0 < (5 < TT cause spontaneous drifts, because accord¬ 
ing to Eqs. 0 and 0 the common frequency of the 
traveling waves obeys 


lim ©1 = lim ©2 = sin 5 

t—^OO t^OO 


K2Gin. ( 10 ) 


Two more equations are obtained from imposing ri 2 = 
0 in 0: 


r2 


-G 2 

Gi cos S 

[{i-4)K2 

ri 

4D 

-Gi 

^ G 2 cos S 

1 

1-1 


( 11 ) 


With ([|) and ( [TT| ) we have three coupled equations 
for three unknowns, ri ^2 and <5. No stationary solution 
with (5 = 7r/2 can be found, but the singularities G 2 ,i = 0 
and ^^ 1,2 = 0 turn out to have a special meaning. In 
particular, if one of the in-coupling strengths A"i ^2 van¬ 
ishes, the corresponding population remains incoherent. 
Numerical continuation around this point shows that in 
order to avoid a negative local order parameter, ri 2 < 0, 
which is unphysical, the whole population transfers to a 
TT-state. In our parametrization this first critical condi¬ 
tion can be written as 


= ±2Ko 


( 12 ) 


One can show that in general the incoherent state, ri 2 = 
0, loses linear stability, if the noise intensity falls below 
a certain value (cf. the Appendix). This happens at 







(AiGAG),, = 8D - 4iGoGo 


(13) 


FIG. 2. (Color online) Bifurcation diagrams for (a) Ko = 
2, Go = 3 and (b) AK = 8, AG = 2. Solid lines from 
Eqs. (12 I, (131, dashed lines with MATCONT [14]. “Zero-lag 
sync” denotes partially synchronous states with zero phase 
lag between the two subpopulations. Parameter regions with 
lag 5 = TT labeled as “yr-states”. Traveling waves (TWl and 
TW2) have non-zero wave speed Q. Insets show enlarged 
areas, (c) Further zooms show coexistence of traveling waves 
with other collective states. Red dotted lines obtained with 
MATCONT, plotted only here. 


Finally, as outlined below, the intersection given by 
G 2 ,i = 0 and (Hi coincides with the origin of bistabil¬ 
ity. Note that the aforementioned conditions are exact. 
Figures l^a) and[^b) depict bifurcation diagrams in the 
planes spanned by the coupling mismatches {AK, AG) 
and the average coupling strengths {Kq, Gq), respec¬ 
tively. Solid lines are given by the critical conditions 
(121 and (13 k Dashed lines are obtained on the basis of 
the reduced system ([^ with the help of MATCONT [TT]. We 
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FIG. 3. (Color online) Order parameters r 2 , ri, r and phase lag |5| as a function of average in-coupling strength Kq. (a) Go = 0, 
(b) Go = 2, remaining parameters in both cases: AK = 8, AG = 2. (c) AK = 3, AG = 10, Go = 2 [cf. Eq. Q]. In all panels, 
dots are o bta ined by integrating the full system Q with N = oscillators. All lines are obtained by numerically solving 
Eqs. In order to unfold the hysteresis in (b), besides forward continuation (>) as in the other cases, a backward 

continuation (<l) is performed. Use of colored regions as in Fig. 


detect branch and limit points, since at all lines one eigen¬ 


value vanishes, except at the ones given by (121, because 
those do not correspond to real bifurcations, but delin¬ 
eate two analogous partially synchronous states: zero-lag 
and TT-states. We additionally test all these findings by 
numerically calculating the eigenvalues of the Jacobian 
of ^ with the fixed points given by ([^-([TT|. When the 
two lines given by (12) and (13) intersect, the boundaries 
(El cease to exist. We emphasize two distinct routes to 
TW states; TWl is surrounded by 7r-states, TW2 by clas¬ 
sical zero-lag sync states. Delimiting lines approach each 
other, see insets for enlarged areas. We further find bista¬ 
bility between incoherence and zero-lag or 7r-states, see 
panels (a) and (b) in Fig. respectively. The bistable 
areas are circumscribed by two lines that intersect at the 
points given by AG = ±2Go and (13l. The location of 
this intersection determines the type of bistability. In¬ 
terestingly, traveling and non-traveling wave states can 
coexist in small parameter regions. We show this in Fig. 
I^c). In particular it is observed that traveling waves 
can coexist with complete incoherence, as well as with 
TT-states and zero-lag partially synchronous states. 


Figure suggests certain conditions for the various 
collective states. In order to observe 7r-states, mixed 
attractive-repulsive in-couplings have to be included. 
Traveling waves surrounded by 7r-states are possible, 
if one includes mixed in-couplings without mixed out- 
couplings. In contrast, traveling waves surrounded by 
zero-lag sync states exist, if there are mixed out-couplings 
and a non-zero mismatch without mixing in the in¬ 
couplings. Bistability between incoherence and zero-lag 
sync requires lack of mixed in-couplings, but the pres¬ 
ence of mixed out-couplings. Finally, bistability between 


incoherence and 7r-states is possible by combining mixed 
in-couplings with vanishing mixing in the out-couplings. 
These conditions appear to supplement consistently the 
observations in Refs. mm- In particular, we verified nu¬ 
merically that the traveling waves surrounded by zero-lag 
sync can also be observed in the setting studied in Ref. 
[^, if a small mismatch in the in-couplings is present, as 
suggested by our bifurcation diagram in Fig. ia). 

It is worth asking how crucial asymmetric interactions 
are for the discordant synchronization patterns discussed 
here. As it is easily seen from Eq. (|^, the interactions 
are symmetric if in- and out-coupling strengths balance 
each other such that the equation 0 = KqAG — GqAK 
holds. This condition can be projected onto straight cuts 
through the parameter space. The demarcations would 
not cross the traveling wave areas [see Fig. Ha) and^b)]; 
instead they would divide the bifurcation diagrams into 
parameter regions that contain both routes to traveling 
waves, i.e. TWl and TW2. In other words, asymmetric 
interactions are needed to get traveling waves. Interest¬ 
ingly however, those straight cuts would go through the 
TT-state regimes. Thus, asymmetry in the interactions is 
not a necessary ingredient to observe 7r-states. This con¬ 
clusion is not evident from the works presented in mm- 


IV. SIMULATIONS 


Additional insights can be gained by numerically solv¬ 
ing the three coupled equations ([^-(11 1 to get ri _2 and 
6. The only subtlety is that one has to factor in the bi¬ 
furcation values previously obtained in order to correctly 
choose between and (§. Solutions are shown in Fig. 
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FIG. 4. (Color online) The spontaneous drift measured by 
common frequency |{0)t|, Eq. ( |10[ | (theory [lines] vs. simu¬ 
lation [dots]) and wave speed |fl|, Eq. ([^ (only simulation), 
(a) Go = 2, (b) Kq = 6; remaining parameters: AK = 8, 
AG = 2 [cf Eq. I^]. Simulation of A = 10'* oscillators, 
averaging over t £ [200, 500], step size dt = 0.01. 


In Fig. I^we compare the common frequency obtained 
from the reduced system, Eq. (10), with numerical sim¬ 
ulations. It serves as an alternative measure to the wave 
speed ^ , which is calculated from the individual instan¬ 
taneous frequencies that do not exist in the analytical 
treatment. Therefore for the wave speed no compari¬ 
son with theory is being made. In Fig. zero-lag syn¬ 
chronous and TT-states become unstable in the regime of 
traveling waves. One can observe that both measures, the 
common frequency and the wave speed, highlight equally 
well the onset of traveling waves. As mentioned before, 
the waves emerge in frequency pairs, meaning that they 
can travel in both directions, depending on realization of 
random numbers. 

In Fig. results of numerical simulations are super¬ 
imposed on smaller theoretical bifurcation diagrams. In 
Figs. [la)-|5]:c ) local and global order parameters are de¬ 
picted, while in Fig. [^d) the wave speed is plotted. For 
each of the 100 x 100 data points in the Kq x Gq grid, 
the equations of motion (Q are integrated with N = 10* 
oscillators and observables are then averaged over time, 
t G [100,50^. Different initial conditions are chosen in 
Figs. me)-!! g) in order to find bistability in numeri¬ 
cal simulations (50 x 50 data points there). Specifically, 
in Figs. |^e)-[^g) the area circumscribed by dashed and 
solid lines shows 7r-states: the local order parameters ri_2 
attain large values, but due to the anti-phase alignment 
given by the phase lag tt, the total order r is small. In 
Figs. [^a)-|^c) this area is filled with incoherence. Hence, 
numerical simulations agree again very well with theoret¬ 
ical results. We remark that bistability between incoher¬ 
ence and zero-lag synchronous states, as predicted by the 
theory, can be analogously found by varying the initial 
conditions (not shown). 


and compared with the results from numerical simula¬ 
tions. For the latter, initial phases are randomly chosen 
from the uniform distribution [— tt, tt]. For each Kq value 
a long-time average is taken over t G [2500, 5000] with 
integration time step dt = 0.01. Upper panels depict 
the order parameters, while in the lower panels the cor¬ 
responding phase lags are shown. The colored regions 
match those in Fig. [2 and discriminate the different col¬ 
lective states. In Fig. ^b) one can see that at the transi¬ 
tion from incoherence to 7r-state the suborder parameters 
abruptly jump from zero to high values in a hysteretic 
manner. For very long time averages it is expected that 
the hysteresis is washed out due to noise-induced jumps 
between the two stable steady states. We do not report 
this here. In Fig. |^c), around Kq « —3, the abrupt 
change in the phase lag and the non-vanishing order pa¬ 
rameters signal extended stability of traveling waves, as 
discussed for Fig. |^c). Such bistable dynamics appears 
to be a promising topic for future studies. In general, the 
abrupt changes and the local minima in the order param¬ 
eters as a function of the average in-coupling strength Kq, 
as presented in Fig. are of vital interest on their own, 
see e.g. Refs. [T^ and m, respectively. 


V. CONCLUSION 

We explored the rich dynamics that emerge from asym¬ 
metric in- and out-coupling strengths among two mutu¬ 
ally globally coupled oscillator populations. As an illus¬ 
trative example we considered identical noisy Kuramoto 
phase oscillators with non-uniform and mixed attractive- 
repulsive interactions. We observed that the two popula¬ 
tions can partially synchronize in-phase and with a con¬ 
stant phase lag to each other. We referred to the latter 
as “discordant synchronization.” The phase lags induced 
spontaneous drifts. As a result, traveling waves were 
formed in which the whole population oscillated with a 
different frequency than the individual units. However, 
in the state of maximal discordance, where the two par¬ 
tially synchronized populations are anti-aligned to each 
other, the spontaneous drifts disappeared. Correspond¬ 
ingly, we revealed two distinct routes to traveling waves, 
one through diametral two-cluster states, the other one 
through classical one-cluster states. Since the latter are 
ubiquitously investigated in the literature, we expect the 
second route to traveling waves to be more prevalent 
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FIG. 5. (Color online) Simulations (colormaps) vs. theory (lines); panels (a) and (b) show local order parameters r 2 and ri, 
while in (c) the global order r is depicted. Panel (d) depicts the absolute value of the wave speed |fl|. Inset shows a zoom-in. 
Simulation performed with initial phases randomly chosen from uniform distribution [—7r,7r]. In panels (e)-(g) simulations are 
repeated with initial phases (()i(0) = 1 Vi. Remaining parameters are the same as in Fig. [1(b). 


in general. Appropriate experimental setups made up 
of two constituents are conceivable, realized e.g. with 
laser systems metronomes [18] or chemical Belousov- 
Zhabotinsky oscillators m- 

With the help of a Gaussian approximation in the in¬ 
finite system size limit we derived a three-dimensional 
system of coupled ODE’s. This reduced system allowed 
a thorough bifurcation analysis and further analytical 
treatment, in excellent agreement with numerical sim¬ 
ulations of a large but finite number of oscillators. We 
found that physically relevant singularities constitute a 
significant part of the bifurcation scenario. We further 
showed which collective states can coexist. Our results 
help to understand the emergence of discrepancies be¬ 
tween individual and collective rhythms, as is observed 
e.g. in neuronal networks [^. Specifically, if the con¬ 
nection strengths were capable to slowly vary in time. 


one could expect temporal patterns reminiscent of the 
high-frequency rhythmic events observed in hippocam¬ 
pal networks |20j . Attractive (positive) and repulsive 
(negative) couplings are often associated with excita¬ 
tory and inhibitory connections among neurons. This 
is reasonable, since positive couplings tend to increase 
synchrony, which is also the case for excitatory connec¬ 
tions in the brain. In contrast, negative and inhibitory 
connections have in common that they tend to decrease 
synchrony m- Exceptions however exist, see e.g. Ref. 

The individual out-coupling strengths considered 
here are particularly suitable to emulate the role of exci¬ 
tatory or inhibitory neurons. This is pointed out in Ref. 
[3 by referring to Dale’s principle, according to which 
a neuron releases the same set of neurotransmitters at 
all its synapses m- Such a comparison would become 
even more applicable by including an excitation threshold 

























































7 


into the system nunsnn. Whether the combination 
of mixed attractive-repulsive interactions on the level of 
in- and out-coupling strengths is experimentally relevant, 
remains an interesting topic for the future. 

Furthermore, future work should approach real net¬ 
works of networks by considering multiple populations 
[25] . see the Appendix for first steps. It would also be 
interesting to extend the present framework towards 
inertia m and imposed phase shifts EZlllH]. Finally, in 
small oscillator populations additional peculiarities can 
be expected [25] . 


This result follows from adapting the derivation (TU] to 
the present case. Above Dc the whole ensemble is inco¬ 
herent, below Dc the incoherent state loses stability and 
partial synchrony is observed. Now we give an alternative 
derivation on the basis of the Gaussian approximation. It 
recovers the exact condition (A.31 and yields an instruc¬ 
tive generalization of the equations in the main part of 
the manuscript. Inserting the Gaussian approximation 
for a ll phase distributions p{(j),t\K,G) into Eqs. (A.l)- 
(A.21, we obtain a set of coupled differential equations for 
the local mean-field amplitudes and mean phases 
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Appendix: Arbitrary number of populations 


Here we discuss the transition from incoherence to par¬ 
tial synchrony in an arbitrary number of interacting pop¬ 
ulations of arbitrary sizes. To this end, oscillators with 
the same pair of in- and out-coupling strengths are again 
grouped into one population. Then the nonlinear Fokker- 
Planck equation for the one-oscillator probability density 
p{(j),t\K,G) reads 


dpi^,t\K,G) d^p{^,t\K,G) 

dt 902 

- -^[{KR{t) sin [e{t) - 

X p{'P,t\K,G)}] . 


The global mean-field amplitude R{t) and phase Q{t) are 
given by 

R(t)e*®(‘) = {{rK',G'{t) G' (A.2) 


The averages ((...)) = f dK' J dG'... P{K', G') take 
into account all in- and out-coupling strengths, K and 
G, via the corresponding joint probability distribution 
P{K,G). In Eq. (A.2) the averages are further taken 
over the local mean-field variables, = 


d4>’e^‘^'p (0', t\K,G) . By a linear stability analysis of 
the incoherent state, p{cj),t\K,G) = l/(27r) VAT, G,t, one 
can exactly derive the critical noise intensity for the syn¬ 
chronization transition: 


{{K'G')) 


(A.3) 


f'iC.G 


- rK,GD + 



K 


Qk,G = 


X {{rK',G'G' cos {eK',G' -Qk,g))), 

^K]G + ^K,G r. 

2 


X {{rK',G'G'sin {Qk\G' -^k,g)))- 


(A.4) 


Let us consider small perturbations drK,G(t) of the in¬ 
coherent state, rK,G{t) = 0 '^K,G,t- The perturba¬ 
tions may give rise to zero-lag synchronous or 7r-states, 
Qk'.G' — ^K,G = WTT, m = 0,1. Accordingly, we can 
separate the network of networks into two groups, 1 and 
2, which contain subpopulations with coupling strength 
pairs {K,G) that lead to the same mean phases Qk,g- 
At the same time, the mean phases of the subpopulations 
in the groups 1 and 2 differ by tt. Linearizing around the 
perturbations frjc .G(f) in the two groups separately, we 
obtain from Eqs. ^A.4): 


[SrK,G] 1 = - [^rK,G]i ^ 

X [{{^'rK-,G-G'))^ - {{5rK-,G'G'))^, 

[5rK,G] 2 = - [5rK,G\2 ^ ~ ^^2 

X [{{SrK',G'G'))^ - {{SrK',G'G'))^]. 

(A.5) 

Here, the two separated networks of networks are la¬ 
beled by the indices 1 and 2. According to Eq. 
the two perturbations can be put together as SR{t)C = 
{{Stk',G'G')) ^ — ((SrK',G'G'))^, where C is some con¬ 
stant coming from an arbitrary global mean phase. As a 
result we obtain 


S'R{t) 


-D+^{{K'G'}} 


6R{t), 


(A.6) 


which leads to the critical condition (A.3l. 


D, 


2 
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